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1. Introduction 

The systematic study of sub-dominant terms in scaling provides researchers with 
guidelines on how to extrapolate e.g. finite-size results to the thermodynamic 
limit [Q, |], ||, ^, ||. Also, it may in itself bring out connections to underlying 
physical properties. These latter may be universal, such as the relationship 
between critical free-energy finite-width correction and conformal anomaly of the 
corresponding universality class [[}| , or otherwise relate to details of the system under 
consideration In the present work we investigate the possible existence of a link 
between lattice structure and presence (or absence) of specific subdominant terms 
in finite-size scaling. A numerical analysis is made of correlation-length and free 
energy data at criticality, as given by the largest eigenvalues of the transfer matrix 
(TM) ||, for Ising and three-state Potts ferromagnets on strips of square, triangular 
and honeycomb lattices with homogeneous, isotropic nearest-neighbour couplings, and 
periodic boundary conditions (PBC) across. 

We recall the following results from conformal invariance J9| for the critical spin- 
spin correlation length, £jv, and (dimensionless) negative free energy per site, /n, on 
strips of width N sites with PBC @, @, [h], |ll): 

A/tt&v = V + «€ + b^N- Ul + . . . (1) 



A 2 (/at - /«,) = ctt/6 + a f N-" + b f N^ +... (2) 

where rj is the decay-of-correlations exponent and c the conformal anomaly number, 
with respective exact values rj = 1/4, c = 1/2 (Ising); r\ = 4/15, c = 4/5 (three-state 
Potts) |(| ^, In terms of the two largest eigenvalues A^ and A^ of the column- 
to-column TM, Cjv 1 = CM A w/ A Jv); N fN = ClnA^; the factor ( is unity for the 
square lattice and, in triangular or honeycomb geometries, corrects for the fact that 
the physical length added upon each application of the TM differs from one lattice 
spacing 0,0, 

Studies of the operator content of conformally invariant theories, and the 
perturbation theory (for finite-size systems, with A -1 as the perturbation variable) for 
the corresponding operator product expansion []l4[ have shown how, in two dimensions, 
the allowed values of the exponents u>, oj\ etc are related to the scaling dimensions 
of the respective set of (model-specific) irrelevant operators. These exponents are 
therefore universal, as the set of their allowed values is fixed for a given model. It 
must be noted that terms with lo = 2, the so-called "analytical" corrections, are 
expected in any theory, as they are related to the conformal block of the identity 
operator |Q. It has also been shown jl5], [l6| that no first-order corrections are 
expected for ground-state energies, the dominant terms being second-order; this will 
be particularly relevant for the three-state Potts model, as seen below. 

Many numerical studies of corrections to scaling pertaining to two-dimensional 
classical (i.e. Ising, Potts etc) systems have actually been carried out in their one- 
dimensional quantum counterparts |l6| , |l7 , 18, 19| ], by taking advantage of well-known 



correspondences [£0|, |21| . While this is expected to have no effect on the determination 
of the universal correction exponents, the amplitudes a^, a/, 6j, 6/ etc are generally 
believed to be non-universal, though e. g. in certain quantum chain systems they 
display similar dependences (Privman-Fisher universality) pT| on the anisotropy 



parameter 1 17|. On the other hand, the relationship between (l + l)-dimensional chains 
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and two-dimensional spin systems is such that the latter are necessarily located on 
a square lattice. Thus, much less is known about corrections to scaling for spins on, 
say, triangular or honeycomb symmetries than on their square counterpart. So far, 
the only explicit (published) reference to the connection between lattice symmetry 
and corrections to scaling seems to be a remark on the fact that, on a square lattice 
operators of spin ±4 will appear, giving rise to iV~ 2 corrections [ H . 

Our main purpose here is to estimate, through numerical work, the first (i.e. 
lowest-order) few amplitudes and exponents as given in Equations (|]J) and (|2|), for 
Ising and three-state Potts spins on square, triangular and honeycomb geometries. We 
shall usually assume, for each model, definite values for the two or three exponents 
used in our fits, taking our hints from conformal invariance theory [fL4[ p"5[ fL6| . 

In particular, we shall seek instances in which, for a given spin model and exponent, 
the corresponding amplitude appears to vanish for one or more lattices and is non- 
zero otherwise. The idea of linking accidentally (or otherwise) vanishing amplitudes to 
underlying physical properties has been exploited fruitfully in the past Q| : the absence 
of vacancy corrections (with an exponent 4/3) has been demonstrated for some Ising- 
like models, implying completeness of the corresponding set of irrelevant operators; 
investigation also showed that this is particular to Ising systems, and that for q = 2 + e 
Potts spins such corrections are of order e. In the present case, we expect to probe 
the interplay between lattice symmetries and the set of irrelevant operators for the 
respective spin systems. To our knowledge, no attempt to draw such a connection 
has been made, despite the wealth of data available for two-dimensional conformally 
invariant systems. 

2. Ising ferro magnets 

We start by considering the Ising model. Exact expressions for the eigenvalues of 
the TM are available for all lattices concerned: square [[j} |2^, £23], triangular pif 
and honeycomb jfl| |2q] . From their inspection, one readily sees that the finite- 
size estimates and Jn of Equations (|l]) and (||) must have well-defined parity 
as functions of TV -1 . Therefore we assume, consistently with: (i) analytical evidence 
derived for £jv orL square ffl and honeycomb |nj lattices (ii) numerical analysis of 
/jv on the square lattice 0]; (iii) the universal result connecting /jv and conformal 
anomaly and (iv) results for quantum Ising chains [|l7|, [l9|], that for all cases one 
has: 

N/n£ N = 1/4 + a^N- 2 + b^N' 4 + ... (3) 

N 2 (f N - foo) = tt/12 + a f N- 2 + 6/AT 4 + ... (4) 

The bulk free energies foo can be calculated in closed form [ |26[ | , again for all three 
lattices. Truncating the series above at 7V~ 4 , one gets finite-size approximants (also 
referred to as two-point fits) to the amplitudes a^j, b^j from pairs of £jv and /at for 
consecutive widths N — 1 and N (N — 2 and N for the honeycomb) . The approximants 
themselves still exhibit a weak iV-dependence, on account of the truncation just 
mentioned. In fact they usually converge rather smoothly as N increases, allowing 
reliable extrapolations to be produced; for example, for the square lattice agrees 
with the exact result Q to two parts in 10 5 . We have found the sequences of two-point 
fits for a£, to behave better than those for a/, bf. For these latter on triangular and 
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honeycomb geometries, as well as on the triangular lattice, numerical instabilities 
or sudden trend reversals arose for large N; the worst such cases were a,f, bf on the 
triangular lattice where behaviour changed abruptly for N > 11. We did not pursue 
the analysis of such deviations, since by then we already had a fairly large sample of 
well-behaved data from which to extrapolate (albeit with less accuracy than in other 
instances, where monotonic trends seem to extend all the way as N — > oo; for instance, 
from a logarithmic plot one finds for the honeycomb a^(N) cx N~ 6A for 6 < N < 46 ). 
Graphical illustrations and extrapolated numerical estimates are shown respectively 
in Figure |l| (where use of 1/N on the horizontal axis is for ease of representation, not 
implying assumption of a specific scaling form), and Table 1. 



able 1. Extrapolated (N — > oo) amplitudes for Ising model (see Equations (|3|) and 
)). Uncertainties in last quoted digits shown in parentheses. 





Lattice 




k 


a f 


b f 


Square 
Triangular 
Honeycomb 


0.102810(2)(t) 

< io- 8 

< IO" 10 


0.2515(5) 
-0.007515(5) 
-0.12022(5) 


0.150730(2) 

< 10" 6 

< 10"° 


0.385(1) 
-0.01165(5) 
-0.1865(10) 



(t) Exact : tt 2 /96 = 0.102808379 ... jlj 



The above results strongly suggest that the coefficients of the N 2 corrections are 
exactly zero in both Equations (^) and (^), for triangular and honeycomb geometries. 



3. Three-state Potts ferromagnets 

In order to check whether the (apparently exact) vanishing of the analytical corrections 
is a model- independent, purely lattice-related phenomenon, we proceeded to study 
the next simplest spin system, the three-state Potts ferromagnet. Exact critical 
temperatures and exponents, and closed forms for bulk free energies, are again 
available for all three lattices |2(jj| . Compared to the Ising case, the main differences 
are: (i) no exact expressions are forthcoming for the eigenvalues of the TM, implying 
that one must rely on numerical diagonalisation, and also that the simple argument 
for definite parity of eigenvalues as functions of TV -1 , plausible in the Ising case, need 
not apply here; and (ii) (corroborating the point just made) a correction exponent 
ujq = 4/5 is expected to arise p7f , overshadowing higher-order terms. 

We tackled (i) by generating the data displayed in Table 2, for strip widths 3 — 14 
(square); 3—12 (triangular) and 4—14 (honeycomb; N even). The data in part (b) 
are displayed without the geometric factor £, in order to match the bulk free energies 
as given in closed form in Ref. (26). Thus they must be multiplied respectively by 
2/V3 (triangular) and V% (honeycomb) to fit Equation (||) with c = 4/5. 

As regards (ii), we took recourse to predictions from conformal invariance, namely: 
(a) analytical (^V -2 ) corrections are always expected [[u]; (b) first-order corrections 
to the ground-state free energy must be absent jl5|, [l(|. Therefore we assumed: 

N/tt£ n = A/l5+a^"N- 4/5 +af J "N- 8/5 +a i N- 2 +af°N- 12/5 + . .. ; (5) 

N 2 (f N - foe) = 2tt/15 + af°N~ 8 ^ 5 + a f N- 2 + A^ 12 / 5 + . . . . (6) 

Because of the small value of loq, second-, and even third-order perturbation 
terms may give rise to corrections of comparable magnitude to the analytical ones. 
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Figure 1. Two-point fits of coefficients in Equations and (|4|). Splines are guides 
to the eye. Lattices and JV-ranges displayed are (top to bottom) square: a^, b^: 
9 — 50, a,f, bf. 9 — 40; triangular: a^: 6 — 27; b^: 6 — 22; a,f,bf. 6 — 11; honeycomb: 
a f , 6 f : 8 - 46; a f , bf. 8 - 22 . 
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Table 2. Finite- width data for 3-state Potts ferromagnet 



(a) r) N = N/tt£ n 



N 


Square 


Triangular 


Honeycomb 


3 


0.292 265 133 729 


0.274124 359 026 




1 


0.282 270 742 335 


0.272 011393185 


0.271 730 743 990 


5 


0.277 362 861763 


0.270 966103 989 




6 


0.274 717 716167 


0.270 319 649167 


0.270 584 629 588 


7 


0.273119 422 538 


0.269 867269140 




8 


0.272 058 600 969 


0.269 527830 301 


0.269 822 520 229 


9 


0.271 304 554 006 


0.269 261539 672 




10 


0.270 741 340 874 


0.269 045 974 456 


0.269 315 733 294 


11 


0.270 304 823 980 


0.268 867303 915 




12 


0.269 956 670152 


0.268 716 440 914 


0.268 957 780 730 


13 


0.269 672 553 348 






14 


0.269 436 295 211 




0.268 691595 289 


oott) 


4/15 


4/15 


4/15 




(b) (l/JV)mA° 




N 


Square'*) 


Triangular 


Honeycomb 


3 


2.121091261980 


2.002 959 900 939 




1 


2.097 704 520 030 


1.985 008 921770 


2.287 656 371786 


5 


2.087460 663 806 


1.976 777367483 




6 


2.082 063 689 903 


1.972 321291068 


2.279 250 879 301 


7 


2.078 863 356 335 


1.969 638 661220 




8 


2.076 806 333 203 


1.967899 052 250 


2.276 305 887939 


9 


2.075 404 683 690 


1.966 707034 799 




10 


2.074 406 246134 


1.965 854 709111 


2.274 943 228 054 


11 


2.073 669 695186 


1.965 224 253 551 




12 


2.073110 712 019 


1.964 744 836 946 


2.274 203 278 951 


13 


2.072 676 417341 






14 


2.072 332 267371 




2.273 757 232 635 


oo(**> 


2.070187162 576 


1.962 224155163 


2.272 522 658 739 



Ct) Conformal invariance |9| 

M Data for N = 3 - 11 available in Rcf. g| 

(**) Evaluated in closed form [M 



Obviously, one cannot simultaneously and reliably fit all coefficients displayed in 
Equations (||) and (|^) from the data of Table 2. We thus resorted to selective 
truncations. Recall that our main goal is to check whether the presence or absence of 
analytical corrections follows the same lattice-dependent pattern found above. 

For correlation- length data we started by trying (al) two-point fits for and 

a| w °, assuming = 0; (a2) same for and a^, assuming a| LlJ ° = 0; (bl) three-point 
fits for a^°, ajr and a^; (b2) three-point fits for a^°, a^"° and a^ , assuming = 0. 
Least-squares fits for varying ranges of N were also performed to the forms used in 
(al)-(b2), always giving similar values of \ 2 ~ 10~ 5 — 10~ 6 per degree of freedom. 

For each lattice, all procedures gave reasonably consistent trends for , no doubt 
because of the small value of u>o, and the wide gap separating it from all other 
assumed secondary exponents; on the other hand, since these latter are so close to 
one another, we could generally extract no clear-cut information from their respective 
fitted coefficients: in all cases, any pair of terms (or single term) among a^ u °, a^, a^ 
could do a reasonable job of standing in for effective corrections when the other(s) 
was (were) assumed absent. Further, we noticed that sequences of finite- N estimates 
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0.022 
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0.018 
0.016 
0.014 
0.012 

0.016 
0.015 
0.014 




Figure 2. Finite-./V estimates of in Equation (g) from different truncation 
procedures (see text). Squares: (al); triangles: (a2); stars: (bl); empty circles: (b2); 
full circles: Equation (pt) . Splines are guides to the eye, joining the vertical axis at 
central estimates from Equation ([?]); see Table 3. Top to bottom: square, triangular, 
honeycomb. 



were much more stable than any of the above if we used an ad hoc form inspired in 
the Ising case, namely 

N/ir£ N = 4/15 + a^'N-^ 5 + a^N' 2 + b^N^ . (7) 

Indeed, in Figure || one sees that, although for each lattice all forms (al)-(b2) yield 
estimates of whose extrapolated values might conceivably coincide (as illustrated 
by the splines), data from Equation (f7j) exhibit the smallest residual ./V-dependence 
of all, leading to the least uncertainties upon extrapolation. Such stability shows up 
also in the estimates of of Equation (0) and, to a more limited extent, in those of 
&£. This contrasts with the corresponding sequences of a 2uJ ° etc in procedures (al)- 
(b2), which generally display a much broader variation. See Table 3 for the respective 
extrapolations against 1/N. 

The small absolute values of for triangular and honeycomb lattices remind one 
of the corresponding case for Ising spins. For (a2) on the triangular lattice, and (bl) 
on the honeycomb, error bars actually include the origin. However, turning now to 
the most regular series given by Equation (R), we recall that the respective estimates 
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are reached by crossing the horizontal axis, on extrapolation of rather monotonic 
sequences; so it seems improbable that by e. g. going to larger widths a change of 
trend would occur, causing error bars to be consistent with = in this case. A 
different possibility, connected with the ad hoc character of Equation (Q), is that there 
may be a systematic error implicit in assuming this particular set of scaling corrections. 
One might ask what combination, if any, of power-law corrections would be consistent 
with the absence of analytical terms. However, having already examined all reasonable 
forms suggested by theory (namely procedures (al)-(b2) and associated least-squares 
fits), and having found none among them associated with clearly superior results, as 
regards smooth convergence and/or quality of fit, we decided not to proceed along 
these lines. 



Table 3. Extrapolated (N —* oo) amplitudes for Potts correlation lengths. See 
text for definitions of procedures (al)— (b2), and Equation Uncertainties in last 
quoted digits shown in parentheses. 



Procedure 


V 


2<jJn 




3^o 


h 






Square 






(al) 


0.0189(7) 


0.0352(6) 








(a2) 


0.0180(5) 




0.119(6) 






(bl) 


0.015(2) 


0.04(3) 


0.05(5) 






(b2) n 


0.016(1) 


0.038(15) 




0.13(5) 




Eq. <|VJ) 


0.01727(4) 




0.133(1) 




0.10(2) 






Trian 


gular 






(al) 


0.0150(2) 


0.000(1) 








(a2) 


0.0149(1) 




0.001(1) 






(bl) 


0.0141(6) 


0.012(7) 


-0.017(10) 






(b2) n 


0.0144(3) 


0.006(3) 




0.012(6) 




Eq. <|2P 


0.01474(1) 




0.0046(1) 




-0.026(10) 






Honeycomb 






(al) 


0.0161(7) 


0.0045(30) 








(a2) 


0.0163(4) 




0.008(3) 






(bl) 


0.0176(20) 


-0.016(25) 


0.028(31) 








0.017(1) 


-0.005(9) 




0.016(14) 




Eq. <|2P 


0.01653(2) 




0.0061(1) 




-0.31(2) 



We now turn to free energies. Attempts to include a first-order term, a^°iV~ Wo in 
Equation (g) consistently produced very small, and steadily decreasing with increasing 
AT, values of (< 10~ 5 , on extrapolation) for all three lattices. Thus, our numerical 
data are in entire accord with the prediction from conformal invariance Jl5| , |l6f that 

Going to higher-order terms, we first recall earlier results, all for the square 
lattice or, equivalently, quantum chains. A single-correction estimate gave u> = 2 Q. 
Assuming terms proportional to N 2uj ° to be negligible, contradictory results were 
obtained upon comparing numerical estimates and conformal-invariance predictions 
for gap amplitudes related to N~ 2 Q. This was later explained by showing that, 
though the amplitude of the TV 2 " terms was indeed small, it could not be neglected. 
Specifically as regards free energies, the following values were calculated, using the 
notation of our Equation (|): a 2 p = 0.001721(3), a f = 0.0280(5) @. 

Again, we tried selective truncations of Equation (^|), namely (al) two-point fits 
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for and aj, assuming ay' = 0; (a2) same for OT and a^ , assuming aj = 0; 
(b) three-point fits for a 2uJ0 , a/ and a^ . Least-squares fits for varying ranges of 
N were also performed to the forms used in (al)-(b). The results were generally 
undistinguished, much as was the case for 77 above. The only exception was (a2) for 
the square lattice where, contrary to all other cases, we found a steadily growing 
with increasing N. This clearly signals, consistently with earlier results [^[ that 
for the square lattice analytical terms are present with a large coefficient. 

In order to make contact with conformal invariance work, we investigated the 
convergence of a 2uJ ° and aj in procedure (al), for the square lattice. We found these 
estimates to vary significantly with N; in fact, Bulirsch-Stoer p8| , |29[ ] extrapolations 
pointed to effective corrections ~ N~ x , x ~ 2.5, with respective final estimates 
a 2u} ° ~ 0.02, a/ ~ 0.25. This gives a value ~ 12.5 for the ratio af/a 2u, ° , to be compared 
to Reinicke's estimate, 16.3(3) ]19| ]. Given the number and severity of approximations 
involved in our calculations, this result may be regarded as broadly consistent with 
universality of the ratio a / /a 2uJ ° , as expected from conformal invariance. 

Since, for all procedures and lattices, the overall picture was very similar to that 
described in the preceding paragraph, we decided to try the simpler scheme of Ref. Q , 
namely assuming effective one-power corrections N 2 (Jn — foo) — 27r/15 + afj N~ u ° , 
and allowing ujq to vary, searching for good fits. We determined the following 
optimum values for ujq: 1.9 — 2.0, (square lattice, again in accord with the early 
estimate iv = 2 |Q); 1.6 — 1.7 (triangular); ~ 1.5 (honeycomb). Similarly scattered 
results arose recently from Monte-Carlo calculations of magnetisation, susceptibility 
and specific heat for three-state Potts spins on square N x N systems J30|, where 
widely differing values of ojq were obtained for the corrections of each quantity. As 
regards the possible connection between lattice symmetry and presence, or absence, 
of analytical corrections, the above one-power corrections are broadly consistent 
with a similar picture to that found for Ising systems. In such scenario, the 
triangular- and honeycomb- effective exponents would reflect a non-zero coefficient 
for N~ 2uJ ° = iV~ 16 , and a null one for the analytical term, while for the square lattice 
the explanation would be as discussed above: both terms are present, with N~ 2 having 
an amplitude one order of magnitude larger than iV -1 ' 6 . 



4. Conclusions 

We have analysed corrections to scaling in critical Ising and three-state Potts 
ferromagnets on the three main two-dimensional lattices. By studying the simplest 
non-trivial systems available and well-defined lattice symmetries, we have left aside 
e.g. anisotropy-induced crossover phenomena and dealt with power-law finite-size 
corrections alone, avoiding marginal operators (which arise in four-state Potts systems, 
for instance) with their associated logarithmic terms. We have given strong evidence 
that the amplitudes of the N~ 2 corrections of inverse correlation length and free energy 
of Ising systems vanish both on triangular and honeycomb geometries, but not on the 
square lattice. 

For the inverse correlation length of the Potts model, the main correction is well 
fitted by an TV -4 / 5 term, in accord with theory and earlier numerical work |l^, |l8]. |i"9[ ]. 
We have found the values of the iV -2 coefficients for the triangular and honeycomb 
lattices to be certainly much smaller than for a square geometry, though we have not 
been able to ascertain that they vanish. 
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For Potts free energies, the N~ 4 / 5 term is always absent, as predicted by conformal 
invariance |L5| , fL6| ; we have estimated the power arising in an effective single-term 
correction to be an apparently lattice-dependent exponent in the range 1.5 — 2.0. Our 
results for the square lattice are consistent with earlier work, in predicting iV -16 and 
TV -2 corrections whose amplitude ratio is ~ 12.5, to be compared with the previously 
obtained estimate 16.3(3) As regards triangular and honeycomb geometries, 

the respective exponents for the effective single-term corrections are much closer to 
2luq = 8/5 (the second-order perturbation value predicted by conformal invariance) 
than to 2. This may, or may not, mean that the N~ 2 term is absent in these cases. 

Several questions arise: (i) can one prove, e.g. on the basis of conformal invariance 
properties, that the amplitudes of the N~ 2 corrections are exactly zero for triangular 
and honeycomb Ising systems? If so, (ii) how do such amplitudes behave e.g. in 
an anisotropic triangular lattice as one crosses over towards a square symmetry, by 
reducing the strength of the bonds along one of its main directions? And (iii) how 
does the 2 + e-state Potts model behave in triangular or honeycomb geometries, as 
regards the N~ 2 correction? (iv) Do the amplitudes of the 7V~ 2 corrections really 
vanish for the three-state Potts model on the same lattices as it appears to be the 
case for Ising systems? If so, can it be established by conformal invariance theory? 
We plan to investigate numerical aspects of these and related matters in future work. 
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